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1. Introduction 

In Ref. p] with S.-H. Tsai, exact results were given for the partition function of the zero- 
temperature g-state Potts antiferromagnet, or equivalently, the chromatic polynomial, 
on open and cyclic chain graphs composed of m repetitions of p-sided polygons connected 
to each other by line segments. In this paper we generalize this work and give the full 
Potts model partition function for these families of graphs. The g-state Potts model has 
long been of interest in the study of phase transitions and critical phenomena [2j. On 
a lattice, or, more generally, on a graph G, at temperature T = l/(fcg/3), the partition 
function for this model is Z = 1 e~ l3n : with the Hamiltonian Ti = — JY"! 
where J is the spin-spin interaction constant, % and j denote vertices on G, is the 
edge connecting them, and <7j are classical spins taking on values in the set {l,...,q}. 
We use the notation K = (3 J and v — e K — 1. Thus, for the Potts ferromagnet (J > 0) 
and antiferromagnet (J < 0), the physical ranges of v are v > and — 1 < v < 0, 
respectively. For the Potts antiferromagnet (PAF), J < so that, as T — > 0, K — > — oo; 
hence, in this limit, the only contributions to the PAF partition function are from spin 
configurations in which adjacent spins have different values. The resultant T = PAF 
partition function is therefore precisely the chromatic polynomial P(G, q) of the graph 
G, which counts the number of ways of assigning q colors to the vertices of G subject 
to the condition that no two adjacent vertices have the same color. (These are called 
proper g-colorings of G.) 

In general, a graph G = (V,E) is defined by its set of vertices (sites), V, and its 
set of edges (bonds), E. We denote the number of vertices of G as n = n{G) = \V\ and 
the number of edges of G as e(G) = \E\. The families of graphs to be considered here 
are open and cyclic chains of polygons connected by line segments. One may regard the 
chain as being oriented so that the longitudinal direction is horizontal. Each polygon is 
connected to the chain at two vertices, such that there are e\ edges of the polygon above 
the chain and e 2 edges below the chain, and there are e g edges between each polygon 
(where g stands for gap). Some illustrative examples are given in Fig. [T](from [1]). The 
basic subgraph unit of the chain is thus a polygon with 

P = e 1 + e 2 (1.1) 

edges, connected to a line segment with e g edges. The full chain with open (o) 
or cyclic (c) boundary conditions (BC) is comprised of m repetitions of this basic 
subgraph comprised of the p-gon and e 9 -length line segment, and is denoted, as in 
PQ, by G , ei,e 2 ,e 9 ,m;BC- Since the two sides of the chain are equivalent, the interchange 
ei -H- e 2 leaves it invariant, so G ei ^ 2>eg>m -Bc — G e2 ^ ei>eg ^Bc- This implies that all of the 
quantities to be presented below are also invariant under this interchange. Indeed, in 
some of these quantities, the numbers e\ and e 2 only enter in the form of their sum, p. 
Since the (m + l)'th member of a family of strip graphs of this sort can be obtained from 
the m'th member by gluing on an additional basic subgraph unit or, in the case of the 
cyclic strip, by cutting the strip transversely, inserting an additional basic subgraph unit 
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and regluing, these are recursive families, in the sense of [3]. The numbers of vertices 
and edges of these graphs are 

n(Gei,e2,e g) m; ) = (P + Zg ~ + 1 , (l- 2 ) 

w(G ei ,e 2 ,e s ,m;c) = (P + e g ~ l ) m , (l- 3 ) 

and 

e ( G 'e 1 ,e 2 , e9 ,m;o) = e(G eue2 , e ^ m . c ) = (p + e 9 )m . (1.4) 

Clearly, G eue2>eg , m . i0 and G , ei ,e 2 ,e 9 ,m;c are planar graphs. In the context of statistical 
mechanics, one takes e\ > 1 and e2 > 1, whence p > 2, since setting either ei = or 
e2 = would mean that a spin <jj would interact with itself rather than with neighboring 
spins. However, in the context of mathematical graph theory, one may formally consider 
the case where e\ — 0, — 1 or e\ = 1, e2 = 0, whence, p — 1. In these cases, the 
graphs G ei;B2te m; BC contain loops, where a loop is defined as an edge that connects a 
vertex back to itself. 

One motivation for the present work is to understand how the results of [T] can be 
generalized to finite temperature and ferromagnetic as well as antiferromagnetic spin- 
spin couplings. Another is to get further insight into how properties of a graph affect 
the Potts partition function Z or equivalent nTutte polynomial (see further below). A 
particular appeal of the chain graphs considered here is that the results are sufficiently 
simple that one can study them in considerable explicit detail. For recursive strip 
graphs of length m basic subunits, Z is a sum of m'th powers of certain algebraic 
functions, generically denoted as A's. Although calculations of Z have been done on 
wider strips of regular lattices, as the strip width increases, the results rapidly become 
quite complicated. For example, for the cyclic (or Mobius) strip graph of the square 
lattice of width L y = 3 and length L x = m, there are Q = 20 different A's, and many of 
these are solutions of algebraic equations of sufficiently high degree so that they cannot 
be expressed in closed analytic form plj. Although for a given width, fewer A's occur 
for an open strip than for a cyclic strip, even for width L y = 3, there are five of these, 
including four that are solutions of a 4'th order algebraic equation, rendering an explicit 
expression rather cumbersome [5j |6]. Thus, it is valuable to investigate the effects of 
graphical properties on Z for families of strip graphs where one can obtain explicit exact 
closed-form analytic solutions for the A's that enter. 

2. General Background 

In this section we discuss some general background material relevant to our study. Let 
G' = (V, E') be a spanning subgraph of G, i.e. a subgraph having the same vertex set 
V and an edge set E' C E. Then Z can be written as [7] 

Z(G,q,v)= q k[G,) v e{G,) , (2.1) 
G'CG 
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(a) 




Figure 1. Illustrations of cyclic and open polygon chain graphs G eil e 2 ,e ,m and 
G ei> e 2 ,e 9 ,m ;o with(ei,e 2 ,e s ,m) =(a) (2,2,0,6), (b) (2,2,1,4), (c) (2,3,2,3). For the cyclic 
(open) chain graphs, the rightmost vertex on each graph is identified with (distinct 
from) the leftmost vertex at the same level, respectively. 



where k(G') denotes the number of connected components of G' . As is evident from 
( 12. ip . Z(G,q,v) is a polynomial in q and v. In the ferromagnetic case with v > 0, 
Eq. (12. ip allows one to extend the definition of q from the positive integers to the 
positive real numbers while maintaining Z(G,q,v) > and hence a Gibbs measure. 
Since k{G') > 1 V G', it follows that Z(G,q,v) always has a factor of q. It is thus 
convenient to define the reduced partition function Z r (G,q,v) as 

Z(G,q,v) = qZ r (G,q,v) (2.2) 

For a graph G, let us denote G — e as the graph obtained by deleting the edge e and G/e 
as the graph obtained by deleting the edge e and identifying the two vertices that were 
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connected by this edge of G. This operation is called a contraction of G on e. From Eq. 
(12. ip . it follows that Z(G,q,v) satisfies the deletion-contraction relation 

Z(G, q, v) = Z(G - e, q, v) + vZ(G/e, q, v) . (2.3) 

The Potts model partition function is equivalent to the Tutte polynomial, 
T(G, x, y), an object of considerable interest in mathematical graph theory. For a graph 

G [5115], 

T(G, x,y)= ( x ~ l) k(G ' yk[G) (y - 1) C(G ' } , (2.4) 

G'CG 

where c(G') denotes the number of (linearly independent) cycles in G' . Note that 
c(G) = e(G) + k{G) — n(G). As is clear from (12. 4D . T(G,x,y) is a polynomial in x 
and y. All of the families of graphs considered here are connected, so that k(G) = 1. 
Let us define 

x = l + -, y = v + l, (2.5) 

v 

so that q = (x — l)(y — 1). Then the equivalence between Z(G,q,v) and T(G,x,y) is 
given by 

Z(G,q,v) = (x- l) k ^ G \y - l) n ^T(G,x,y) . (2.6) 

The special case v = — 1 defines the T = Potts antiferromagnet, and sets the variables 
in the Tutte polynomial equal to x = 1 — q and y = 0; in this case these functions yield 
the chromatic polynomial: 

Z(G, q, -1) = P(G, q) = {-q) k ^ (-l) n ^T(G, l-q,0). (2.7) 

If G contains a loop, then T(G, x, y) contains a factor of y and P(G, q) vanishes 
identically, since it is not possible to satisfy the proper g-coloring condition. 
We recall two elementary results. For any tree graph T n with n vertices, 

T(T n , x, y) = x n ~\ Z(T n , q, v) = q(q + v) n ~ l . (2.8) 

For the circuit graph with n vertices, C n , 

T(C n ,x,y) = ?-^ , Z(C n ,q,v) = (q + vr + (q-l)v n . (2.9) 
x — 1 

where c\ = xy — x — y — q — 1. Note that although T(C n , x, y) is expressed in Eq. (12. 9p as 
a rational function, it is actually a polynomial, as is guaranteed by Eq. (12.41) . Explicitly, 
T(Ci,x,y) = y and, for n > 2, T(C n ,x,y) = y + ^j=i xJ - A similar comment applies 
to expressions below involving (x — 1) denominators. For the special case v = —1, since 
P(C n ,q) has a factor q(q — 1) for n > 2, it is convenient to define the polynomial D n 
via 

»-f^g<-i)f:y 2 -*- (2 - io) 

(For n — 1, P(Ci, q) = 0, since C\ is a vertex with a loop; hence also D\ = 0.) 
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3. Calculations and Results 



Using a systematic application of the deletion-contraction theorem, we have calculated 
the Potts partition function and equivalent Tutte polynomial for the open and cyclic 
chain graphs G eue3)egtTn . and G eue3)egtTn;c . For the open chain graph G eue ^ egim]0 we find, 
for the Tutte polynomial, 

T{G eue2>egim . !O ,x,y) = (\ T;0 ) m , (3.1) 

where 

At,o = T(T eg+1 , x, y) T(C P , x, y) = x^ (^p) ■ (3.2) 

Note that At,o only depends on e% and e 2 via their sum, p. In general, if a graph G 
can be expressed as G\ U G2 such that G\ fl G2 is a single vertex, then T(G,x,y) = 
T(G\,x,y)T(G2,x,y). This factorization property holds for the basic subgraph units 
of G eije2iegjm;0 , since the p-gon intersects the line segment at a single vertex. The 
factorization property also holds, a fortiori, for the open chain G eiie2 , egim;0 itself. These 
facts imply the properties that (i) At,o has the form of a product of the Tutte polynomials 
of the tree graph T £g+ i and of the circuit graph, C p , and (ii) T(G eije2iegim;0 , x, y) has the 
form of a power of a single polynomial, as given in Eq. (13. ip . Corresponding comments 
apply to Z(G eije2iegiTn . ,q,v). 

For the cyclic chain graph G ei ,e 2 ,e 9 ,m;c we calculate 

T(G ei>e2tegtm;c ,x,y) = -l_[(A Ti0 r + Cl (A Til r] , (3.3) 

where At,o and c\ were given above and 

A Tjl = T(C ei ,x, y) + T(C e2 , x, y) + l-y. (3.4) 

It is noteworthy that At,i is independent of e g and only depends on the properties of 
the polygons, as encoded in their edge numbers e x and e 2 . 

From these general results, one can consider special cases of the various graphical 
edge numbers e±, e 2 , and e g . For example, for the open and cyclic chain graphs G , 2,2,e !? ,m;o 
and G2 1 2,e g ,m;ci illustrated (for e g = 0, 1 and m = 6,4) in Fig. [TJ, we have 

T(G 2 ,2,e g ,m;oi X, V ) = [x^ (x + X 2 + X 3 + V )] m (3.5) 

and 

T(G 2 ,2,e g , m ;c, x, y) = —- \{x e9 {x + x 2 + x 3 + y)} m 
x — 1 L 

+ (xy - x - y)(2x + y + l) m . (3.6) 

Various evaluations of T(G eije2j£gtm]0 , x, y) and T(G eiie2iesim;c , x, y) for special values 
of the arguments x and y are of interest. In particular, for x = 1, we have 

T(G eue2 , eg , m;0 ,l,y) = (p + y-ir (3.7) 

and 

T(G eue2>eg>m]C , l,y) = (p + y- l)^ 1 (me g +y-l)(p + y-l) 



5 



+ y \P(P ~ !) - e i( e i - !) - e 2(e 2 - 1)) 



(3.8) 



4. Potts Model Partition Functions 



The equivalent Potts model partition functions are 

q,v) = q{\z,o) m 

and 

Z(G eue2 , eg , m;c ,q,v) = (X z , ) m + (q- l)(Xz,i) m , 

where 

A z , = Z r (T eg+1 , q, v) Z r (C p , q, v) = (q + t^g" 1 (q + vf + (q - l)v p 



(4.1) 
(4.2) 



(4.3) 



and 



v e2 Z r (C ei ,q, v) + v e ^Z r (C e2 , q, v) - v* 



v eg 



v e2 \ (q + vf 1 + {q- l)v ei \+v ei \(q + vf 2 + (q - l)v 



,<-'2 



(4.4) 



We note that 



Xz,o = Az,i = v^- 1 (p + v) at q = 



and 



Az i = u 



e„+l 



2(u + l) -v 



p-i 



for ei = e 2 = 1 . 



(4.5) 
(4.6) 



5. Chromatic Polynomials 



For the special case v = —1, our results for Z(G eue2tegtrn;o ,q,v) and Z(G eii e 2jeffim;c , g, i>) 
reduce to to the chromatic polynomials given (for e\ > 1 and e 2 > 1) in [T], 

P(G ei , e2 , eg , m;o ,g)=g(Ap, ) m (5.1) 

and 

P(G eil e 2 , w , q) = (Ap, ) m + (ff - l)(Ap,i) m , (5.2) 

where 

Apo = (g - l) efl+1 A> (5.3) 

and 

A Pil = (-1)^ g- 1 [(1 - g) 6 * + (1 - qf 2 +q-2 
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1) 



p+e a 



''i 



t'2 



(5.4) 



s=2 v ' s=2 

(If either ei = or e 2 = 0, the graphs contain one or more loops and the chromatic 
polynomial vanishes identically.) 



6. Free Energy 

It is of interest to remark on some thermodynamic properties of the g-state Potts model 
on the infinite-length limit of these polygon chain graphs. Although our results for the 
partition function apply for general q, we shall restrict our attention here to integer 
q > 2. We denote the m — > oo limit of the chain graph G e i,e 2 ,e 9 ,m,.BC as {G e l,e 2 ,e 3 ;,Bc} • 
The reduced, dimensionless free energy per vertex for the Potts model in this limit is 
independent of boundary conditions, so we drop the BC subscript. We find 



f({ G e 1 ,e 2 ,e g },q,v) 



lim — In 

m— >-oo fi 



Z(G ei 



e2,e g ,m;BC ) 



(p + e g -l) 
1 

(p + e g -l) 



ln(A 



zo) 



In 



(q + vf* {(q + vf + (q - l)v p } 



(6.1) 



(The actual Gibbs free energy per vertex is Q = —ksTf.) As was discussed in p[], if one 
were to start with a Potts model partition function Z(G,q,v) for a cyclic strip graph 
G with variable real (positive) q, then, for a set of special (s) values of q, denoted {q s }, 
one could encounter the noncommutativity 

> 1/n ^ lim lim Z{G,q,v) 1/n . 



lim lim Z(G, q, v) 

n— >oo q— >q a 



q—tq s n— >oo 



(6.2) 



For G 



ei,e2,e g ,m;c) 



{q s } = {0, 1}. Since we are only interested in (integral) q > 2 here, we 
do not encounter any such noncommutativity. 

Because these chain graphs are quasi-one-dimensional, the Potts model (with either 
sign of J) does not exhibit any critical behavior at nonzero temperature. However, it 
is of interest to investigate how thermodynamic quantities depend on the graphical 
parameters e 1; e 2 , and e g . A first observation is that the reduced free energy / and 
hence quantities that are obtained as derivatives of / with respect to temperature T, 
such as the internal energy per site, U = —df/d/3 and specific heat per site, C = dll / dT, 
only depend on e% and e 2 through their sum, p. This property follows from the fact that 
(in the limit m — > oo) / is determined completely by the dominant A (i.e., the A with 
the largest magnitude), namely \z,o, and this only depends on e\ and e 2 via their sum, 
p. As an explicit example, the internal energy per site is 



U 



J(v + 1) 



(p + e g 



q + v 



+ p 



(q + vy- 1 + (q-l)v p - 1 
(q + v)p +(q- 1)vp 



(6.3) 
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In comparison, we recall that for the Potts model on the line, U\d = —J(v + 1)/ (q + v). 
In the expression (16 .3p for U, one thus sees the interplay of the underlying circuit graph 
including the m line segments, each of e g edges, with the m p-sided polygons, in the 
limit as m — > oo. The limit of Eq. (16. 3p as the temperature T — > oo, i.e., v — > 0, is 

J ( P + e g 



U = ^ 9 at v = . (6.4) 

q \p + e g - I J 

Since (p + e g )/(p + e g — 1) > 1, U(v = 0) is more negative than the analogous infinite- 
temperature limit of U for the Potts model on the line, namely U(0)id = —J/q- One 
can also consider the zero-temperature limit. For the antiferromagnet, this reverts back 
to the analysis of the chromatic polynomial, as in pQ. For the ferromagnet, as j3 — > oo, 

lim U — — J ( P + eg ) for J > . (6.5) 



,s^oo \p + e g — 1, 

Again, this is more negative than the corresponding expression for the Potts ferromagnet 
on the line, which is lim^^oo Uio = — J. In both the high- and low-temperature limits, 
these differences can be attributed to the additional spin-spin interactions due to the 
combination of the repeated polygons attached to the underlying global circuit graph. 
Finally, we note some limiting cases: 

lim /({G eiie2ie }, q, v) = ln(g + v) for fixed finite p (6.6) 

and 

lim f({G e , e „ e }, q, v) = hi(q + v) for fixed finite e„ . (6.7) 

In both of these cases, the dimensionless free energy per site thus reduces to that of the 
n — > 00 limit of the Potts model on a line, / = ln(q + v). This reduction can be ascribed 
to the dominance of one of the two structural parts of the graph - the global circuit in 
Eq. (16. 6p and the polygons in Eq. (\6.7L 



7. Zeros of the Partition Function and Locus B 



It is also of interest to consider the locus of zeros of the partition function. In the limit 
m — > 00, some zeros may merge to form curves denoted as the locus B. It is convenient 
to consider this locus in the q plane for fixed temperature variable v. For the case 
v = —1, i.e., the T = Potts antiferromagnet, B was studied in [TJ. Here we can use our 
general results to analyze the case of the Potts antiferromagnet at temperatures above 
T = and the Potts ferromagnet. Thus, we consider B for v in the full physical range, 
v > -1. 

For the open polygon chain graph, since Z{G ei ^ 2figtm - 01 q, v ) involves just the m'th 
power of a single A term, it has a fixed set of discrete zeros, independent of m, with 
no continuous locus, i.e., B — 0. For the m — > 00 limit of the cyclic polygon chain 
graph, B is a nontrivial locus, determined by the equation \\z,o\ = |Az,i|- For this limit, 
with v = — 1, several properties of B were derived in [I], namely, that B (i) is compact, 
(ii) passes through q = 0, and (iii) encloses regions in the q plane. With the notation 



8 



e s and denoting the smaller and larger of e\ and e 2 , it was shown in pQ that (iv) 
if e s = 1, then $ is the circle \q — 1| = 1, independent of the values of and e g , so 
that q c = 2, where g c is defined as the maximal point at which B crosses the real axis 
(which it always does for the m — > oo limit of these cyclic polygon chain graphs). Two 
additional properties concerning q c values for v — — 1 (denoted (B5) and (B6) in pQ) 
were also shown. 

We use the notation '\ixn m ^. 00 G ex ^ >e m . c = {G ei>e2te . c }. For this m — > oo limit of 
the cyclic polygon chain graph, we find that properties (i)-(iii) continue to hold. The 
proof of the compactness property (i) is a generalization of the proof for v = — 1 given 
in pQ; one uses the fact that a necessary and sufficient condition that B is noncompact 
in the q plane, passing through 1/q = 0, is that the equation \Xz,o\ = \Xz,i\ has a 
solution for 1/q = 0. To show that this is not the case, we extract a factor of q p+e 9~ 1 
from both sides of this equation and define Xz,d = <Z p+es_1 Az )( 2 for d = 0,1. Dividing 
both sides by the factor |gP +e 9 _1 | we have \Xz,o\ = \Xza\- Taking 1/q — > 0, we see 
that this equation cannot be satisfied because Xz,o — > 1 while Xz,i —> 0. This proves 
the compactness of B for general v. Property (ii) follows because, as Eq. (14.51) shows, 
\Xz,o\ = |Az,i| at q = 0. To show property (iii), we evaluate / in Eq. (16.11) (away from 
q s values so as to avoid the noncommutativity). The result is that in region Ri, which 
contains the semi-infinite real intervals q > q c and q < 0, Xz.o has a larger magnitude 
than Xz,i, so / = (p + e g — l)" 1 In Xz,o, as in Eq. (16. lft . while in a region R 2 containing 
an interval of small positive q neighboring the origin, Xz,i has a larger magnitude than 
Xz,o, so / = (p + e g — l)~ 1 lnA^ i i. If these regions were not completely separated by 
the nonanalytic boundary B, then one could analytically continue / from one to the 
other, but this would lead to a contradiction, since / has a different functional form 
in these two regions. This proves that B separates the q plane into regions. Property 
(iv) for v — — 1 was a special consequence of the complete intersection theorem for 
chromatic polynomials, namely that if a graph G is the union G = G± U G 2 such that 
the intersection G\ fl G 2 = K r , where K r is the complete graph on r vertices (i.e, 
the graph in which each vertex is connected to every other vertex by an edge), then 
P(G,q) = P{G x ,q)P{G 2 ,q)/P{K r ,q), where P(K r ,q) = JY^Q ~ s )- This theorem 
does not apply for v ^ — 1, so the property (iv) for v = —1 no longer holds for v > — 1. 

We have obtained a number of results on q c for the general case v > — 1. As we 
have discussed above, for {G eij£2teg - C }, the locus B crosses the real q axis at q = and 
at a maximal point denoted q c . (Depending on the values of ei, e 2 , e g , and v, B may 
cross the q axis at additional points besides q = and q c .) Some exact results on q c are 

q o ({Gi,i, 0iC }) = -2v(v + 2) (7.1) 

q c ({G lxl . c }) = -v(v + 3) (7.2) 

and 

q c ({G 2Xeg ., c }) = -2v Ve r (7.3) 
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For most cases, g c is a root of a higher-order equation. For example, for {Gi 5 i 2;c }, q c is 
the maximal real root of the cubic equation 

q 3 + v(v + 4)g 2 + v 2 (2v + 5)q + 2v 3 {v + 2) = . (7.4) 

These calculations may be compared with our previous results for q c for the infinite- 
length limits of some other strips with periodic longitudinal boundary conditions. 
Denoting lim n _ s . 00 C n as {C}, one has q c ({C}) = —2v, the same as g c ({G2,2,e 9 ; C })- For 
the cyclic or Mobius square-lattice ladder strip pE], 

q c ({sq lad}) = —v(v + 3) , (7.5) 

the same as q c ({Gi,i,i;c})- F° r the cyclic or Mobius ladder strip of the triangular lattice 
[TTj . and also the self-dual (sd) square-lattice ladder strip, 

v(2v + 5) 

q c ({tri lad}) = q c ({sq s d lad}) = — — . (7.6) 

1) \~ — 

These results for q c ({G}) have the general property that, as v increases from —1 to 0, 
i.e., as the temperature increases from to infinity for the Potts antiferromagnet, they 
decrease monotonically from their v — — 1 values to 0, and in the ferromagnetic region, 
v > 0, they are negative. One may also consider the accumulation locus of partition 
function zeros in the v plane for fixed q. However, this is of somewhat less interest than 
B in the q plane for fixed v , since these polygon chain graphs are quasi-one- dimensional 
and hence (for q > 2) B does not cross the real v axis at any point corresponding to 
nonzero temperature. 



8. Flow Polynomials 

An interesting problem in graph theory is the task of enumerating discretized flows on 
the edges of a (connected) G that satisfy flow conservation at each vertex, i.e. for which 
there are no sources or sinks. The flow on each edge can take on any of q values modulo 
q, (so q = Omodg). One arbitrarily chooses a direction for each edge of G and assigns a 
discretized flow value to it. The value zero is excluded, since it is equivalent to the edge 
being absent from G; henceforth, we take a g-flow to mean implicitly a nowhere-zero 
g-flow. The flow or current conservation condition is that the flows into any vertex must 
be equal, mod q, to the flows outward from this vertex. These are called g-flows on G, 
and the number of these is given by the flow polynomial F(G, q). This is a special case 
of the Tutte polynomial for x = and y = 1 — q: 

F(G,q) = (-iy^T(G,0,l-q) . (8.1) 

A bridge on a graph G is defined as an edge with the property that if it is deleted, this 
increases the number of connected components of G by one. If a (connected) graph 
G contains any bridge, then the flow polynomial vanishes identically. The open strip 
graph Gei,e 2 ,eg,m;o contains at least one bridge if e g > 1, and therefore it does not allow 
any g-flows: 

F{G eue2tegtm . , q) = if e 9 > 1 . (8.2) 
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For e g = 0, we find 

F{G eue2 ^ m , , q) = (q- l) m if e g = . (8.3) 

Note that this is equal to [F(C P , q)] m , showing that the flows occur independently 
in each p-gon circuit graph. Algebraically, this follows from the factorized form of 
T(G ei}tlaie m . , x, y) in Eq. (13. ip . This is also clear since, if one considers the flows in the 
p-gons forming the ends of the open chain, these have nowhere else to flow, and, since 
that is the case, the same is true for the flows in all of the p-gons that form the interior 
of the chain. 

Using our calculation of the Tutte polynomial of G ei>e2ie iTn;c , we find, for the cyclic 
strip, 

f(g-l)(g-2) m ife„>l 
F(G eie2e mc ,q) = < 9 ~ . (8.4) 

1 ' 2 ' 9 ' ' ' \(q-l)(q-2) m + (q-l) m if e g = 

Thus, both for e g = and for e g > 1, the cyclic polygon chain graph allows more q 
flows than the open polygon chain graph. This follows because of the freedom of the 
flows in the cyclic case to make a global circuit around the chain. These calculations for 
the cyclic polygon chain may also be compared with the result F(C n , q) = q — 1. One 
sees that if e g = 0, then there are more flows on the G eite2!e m]C graphs, owing to the 
possibility of flows within each polygon. However, if q = 2 and if e g > 1, then, owing to 
the degree-3 vertices where the polygons connect onto the line segments in the chain, 
no g-flows can occur, since the flow conservation condition cannot be satisfied at these 
vertices. 

Since these results hold for arbitrary chain length m, one may consider the limit 
of infinite length and define, as in Ref. [12], a function <fi representing the number of 
g-flows per face of G in this limit, 

<P({G ei ,e 2 ,e g ;BC,q) = 1™ [F(G ei ,e 2 ,e g ,m;BC , (?)] 1/MG) , (8.5) 

m— >oo 

where fc(G) denotes the number of faces of G. Here, fc{G eiye ,^ egtm - ) — m + 1 and 
/c(Gr eije2> e 9im;c ) = m + 2. Unlike the free energy, the function does depend on whether 
one uses open or cyclic boundary conditions, since some flows make global circuits 
around the chain in the cyclic case. For q > 2 (so that nowhere-zero q- flows can occur), 
we calculate 

"(i(v,,,i^/)H , t e9 ~t ( 8 - 6 ) 

— 1 if e g = 



and 



, , - 2 if e„ > 1 

<t>({Ge 1 ,e 2 ,e g ;c},q) = { * -r n " (8 ' 7) 

- 1 if e g — 
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9. Reliability Polynomials 



A communication network, such as the internet, can be represented by a graph, with the 
vertices of the graph representing the nodes of the network and the edges of the graph 
representing the communication links between these nodes. In analyzing the reliability 
of a network, one is interested in the probability that there is a working communications 
route between any node and any other node. This is called the all-terminal reliability 
function. This is commonly modeled by a simplification in which one assumes that each 
node is operating with probability p no de and each link (abbreviated £) is operating with 
probability pg. As probabilities, p no de and pe lie in the interval [0,1]. The dependence 
of the all-terminal reliability function Rtot{G,p no de,Pe) on p no( i e is an overall factor of 
(PnodeT; i-e., Rtot(G, p node , pi) = (p no de) n R(G , p £ ) . The difficult part of the calculation 
of Rtot(G, p no de, Pe) is thus the part that depends on the links, R(G,pg). The function 
R(G,p e ) is given by 

i2(G,ft) = X)p? e) (l-ft)^-^ 

GCG 



(9.1) 



where G is a connected spanning subgraph of G. Each term in this sum is the 
probability that the communication links E e G are functioning (equal to Pg^) times 
the probability that the other links, E—E, are not functioning (equal to (1— pg) e ^~ e ^). 
From its definition, R(G,pg) is clearly a monotonically increasing function of pn G [0, 1] 
with the boundary values R(G, 0) = and R(G, 1) = 1. R(G,pg) is given in terms of the 
Tutte polynomial, evaluated with x = 1 (guaranteeing that G is a connected spanning 
subgraph of G) and y = yg, where 

1 Pe 



Ve = 



I-Pe 



i.e., v t = yt-\ 



I- Pi ' 



by the relation 



R(G,p e ) = p n t -\\ - Pe) e{G)+1 ~ n T(G, 1, - 



Pe 



(9.2) 



(9.3) 



Using our calculation of the Tutte polynomials for G eii e 2 ,e 9 ,m;o and G 



ei,e2,e g ,m 



;c we 



find 



and 



P V e +ea 1 {p( 1 ~Pe) +Pe} 



R\G ei ,e2,eg,m,oi Pe) 

T><r> n ^ — J-P+ e 9- 1 ) m 

^IWi^.eg.mjoWJ — Pe 



(9.4) 



p(l -p e ) +p e 



m—l 



X 



|me g (l -p e ) +p^||p(l - pe) + Pe} 
+ ^(1 - p,f[p{p - 1) - ei ( ei - 1) - e 2 (e 2 - 1)} 



(9.5) 



In general, R(Ge u e 2 ,e g ,m;c,Pe) > R(G eu e 2 ,e a ,m;o,Pe), with equality only at p e = 0, 1. This 
can be understood as a consequence of the fact that with the cyclic boundary condition, 
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there are more possible communication routes linking two nodes than there are with 
the open boundary condition. We observe that for p e G (0,1), R(G eije2tegjm]0 ,pi) and 
R(Ge 1 ,e 2 ,e g ,m;c,Pe) are (i) decreasing functions of m for fixed e±, e 2 , and e g ; (ii) decreasing 
functions of e g for fixed e%, e 2 , and m; and (iii) decreasing functions of t\ for fixed e 2 , e g , 
and m. These properties can be ascribed to the greater probability of communication 
bottlenecks as the respective parameter, m, e g , or e\ increases with the other parameters 
held fixed. We also observe that for fixed p = e\ + e 2 , e g , and m, these reliability 
polynomials increase as |ei — e2 1 decreases. 

As in [13], in the limit of infinite chain length, m — > oo, one may define a function 
r that measures the reliability per node, as 

r ({ G ei,e2,e 9 },^) = hm [R(G ei ,e 2 ,e g ,m;BC , Pi)} V " • (9.6) 



As was discussed in [13] for other strip graphs, this function is independent of the 
longitudinal boundary conditions, so we drop the BC subscript on the left-hand side of 
Eq. (19. 6p . Clearly, for a general {G}, r({G},pi) is an increasing function of pi G [0, 1] 
with the values r({G},0) = and r({G}, 1) = 1. From our exact calculations above, 
we find 



r ({G ei ,e 2 ,e g },Pe)=Pe\p(l-pe)+P^ 3 ■ (9.7) 

As a comparison, for the infinite- length limit of a line graph L n or circuit graph C n , 
r({C},pi) = r({C},pi) = pi = r 1D (pi). Now for p £ G (0,1) (and p > 2), the factor 
[p(l — pe) + p£]V(p+ e 9-i) > i ; and hence r({G ei<e2jeg },pi) > r 1D (pi). Concerning the 
dependence on ei and e 2 (which only enter in the form of their sum, p) and on e g , we 
find from an analysis of the respective partial derivatives, for a fixed pi G (0, 1), that 
r({G eite2teg },pi) is a monotonically decreasing function (i) of p, for fixed e g , and (ii) of 
e g , for fixed p. Furthermore, 

lim r({G eue2te },p e ) = pi (9.8) 

p— >oo 

and 

lim r({G ei>e2ie },p e ) =pi . (9.9) 

e g — >oo 

Thus, in both of these limits, r({G eite2j£g },pi) reduces to r 1D (pi). 



10. Percolation Clusters 



In this section we use our results to calculate a quantity of interest in the area of bond 
percolation. We first briefly mention some necessary background. Consider a connected 
graph G and assume that the vertices are definitely present, but each edge is present only 
with a probability pi G [0, 1]. In the usual statistical mechanical context, one usually 
considers a limit in which the number of vertices n — > oo. An important quantity is the 
average number of connected components (= clusters) in G. For a given G, we denote 
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this average cluster number per vertex as (k)o- This is given by 

[1/n) Eg, k(G')pf\l - Pe )<G)-e(G') 



(k) 



G 



T, & pt G,) ^-Pi) e{G) - e{G,) 



[l /n) E c ,fc(G> 
T v e{G>) 



(10.1) 



where G' is a spanning subgraph of G, as above, and ve was defined in Eq. f 19 . 2 j) . Hence, 
in the n — > oo limit, the average cluster number per vertex, (&)r G i ; is given by 

df({G},q,v) 

0(J 9=1, v=V{ 

Using Eq. fl 1 . 2 j) with Eq. ( 16. ip , we find, for the infinite- length limits of both the open 
and cyclic polygon chains, the average cluster number 



(k) 



{G} 



(10.2) 



<*>{Gn. 



,} 



(p + e g -l) 
1 - Pi 



(p + e g )(l -p e ) + f e - 1 

j=0 



-.p + eg-l 

This may be compared with the result (A;) id = 1 — pi for the infinite line. We have 

p-i 



(10.3) 



(k)iD - (k){G ei , e2 , eg } 



1 - Pe 

p + e g 



- - J=0 

> , (10.4) 

with equality holding only if pi — 1 or p — 1. Thus, if pn < 1 and p > 2, the average 
number of clusters per vertex is greater for the line than for the infinite-length limit of 
the polygon chain graph (with either set of boundary conditions). 



11. Some Graphical Quantities 

Special valuations of the Tutte polynomial of a graph yield various quantities describing 
properties of this graph. In this section we give these. First, we recall some definitions. 
A tree graph is a connected graph with no circuits. A spanning tree of a graph G is a 
spanning subgraph of G that is also a tree. A spanning forest of a graph G is a spanning 
subgraph of G that may consist of more than one connected component but contains no 
circuits. The special valuations of interest here are (i) T(G, 1,1)= N$t{G), the number 
of spanning trees (ST) of G; (ii) T(G,2, 1) = Nsf(G) the number of spanning forests 
(SF) of G; (iii) T(G, 1,2) = Nqssg(G), the number of connected spanning subgraphs 
(CSSG) of G; and (iv) T(G, 2, 2) = N SSG (G) = 2< G \ the number of spanning subgraphs 
(SSG) of G. For both the open and cyclic strips, the last of these quantities is directly 
determined by Eq. ( II. 4p (without the necessity of calculating the Tutte polynomial) to 
be 

NsSG(G ei , e2 , eg , m ;o) = Nss G (G ei , e2 , eg , m . c ) = 2^ m . (11.1) 
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We evaluate our general results for T(G eue2tegtm;o , x, y) and T{G ei ^ egmci x, y) to obtain 
the quantities (i)-(iii). For the numbers of spanning trees, we find 

NsT{G ei ,e2,e g ,m;o) = 

and 

1 



„m— 1 



NsT(G ei: e 2 ,e g ,m;c) = ™p 

The numbers of spanning forests are 

NsF(G ei ,e 2 ,e g ,m;o) = ft* (2> - l) 

and 

V s (2 P - 1) 



P e g+^\pi.P- l )- e ii. e i- l )- e 2{.e 2 -l) 



2 ei + 2 e 2 _ 2 



ei,e2,e g ,m;cj 

The numbers of connected spanning subgraphs are 

NcssG(G ei , e2 , eg , m ; ) = (p + ir 

and 

N C SSG{G ei ,e 2 ,e g ,m;c) = (p + l)™" 1 \(me g + l)(p + 1) 



11.2) 
(11.3) 

11.4) 
11.5) 
11.6) 



rn 



+ y 1 P(P - l ) ~ e i( e i _ !) _ e 2( e 2 - 1) 



;il.7) 



We also give another evaluation of the Tutte polynomial. For a general connected 
graph G = (E, V), one can define an orientation of of G, i.e., a directed graph G = (V, E) 
by assigning a direction to each edge e G E. There are 2 e ^ of these orientations. Among 
these, an acyclic orientation of G is defined as an orientation that does not contain any 
directed cycles. Here, a directed cycle is a cycle in which, as one travels along the 
cycle, all of the oriented edges have the same direction. The number of such acyclic 
orientations is denoted a(G) and is given by the evaluation of the Tutte polynomial 
with x = 2 and y = [14J: 

a(G) = T(G,2,0) . (11.8) 

Equivalently, this is obtained by the evaluation of the chromatic polynomial at q — — 1: 
a(G) = (— l^^P^G, —1). From the results in p], in agreement with our calculations 
here, we have, for the numbers of acyclic orientations of the open and cyclic polygon 
chain graphs 



a (G ? ei,e 2 ,eg,m;o) — 2 9 (2 P 2) 



and 



a(G ei 



e2,e g ,m;cj 



2 e «(2 p -2) 



2 ei + 2 e 2 _ 3 



(11.9) 

11.10) 
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12. Conclusions 

In conclusion, in this paper, generalizing our previous results on chromatic polynomials 
with S.-H. Tsai in p], we have presented exact calculations of the Potts model partition 
functions and equivalent Tutte polynomials for a class of polygon chain graphs with 
open and cyclic boundary conditions. We have evaluated special cases of these 
results to compute the corresponding flow polynomials, reliability polynomials, and 
various quantities of graph-theoretic interest, and have analyzed the dependence on the 
parameters e±, e^-, and e g characterizing the families of graphs. 
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